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.SYMBOL MEANINGS 


E 

Energy 

X 

Position 

Ux 

Cosine of the polar angle with 
respect to the X direction 

t 

Time 

At 

Time interval 

v 

Speed 

S 

Pathlength 

Z 

Atomic number 

C 2 

Speed of light 

m o c 

Rest mass energy 

AS X 

X displacement 

AE 

Energy interval 

P 

Apportioning fraction 

AX 

Position interval 


P 


P(u x |uj.) 


*1 (u) 

G i Cs) 

N 


Particle density 

Cosine of the polar angle of in- 
cidence with respect to the X 
direction 

Angular probability density 
function with respect to u 
given an incident direction u . 
independent of incident and 1 
scattered azimuthal angle 
Legendre Polynomial of index l 
Goudsmit-Saunderson distribution 
expansion terms 
Number of scatterers per unit 
volume 

Polar angle of scattering rela- 
tive to an arbitrary incident 
direction 


♦This work was performed under the 
Commission. 
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SYMBOL MEANINGS 
(continued) 


ct(9 s /S) 

T n (x) 

Au x 

J(X,E, Ux ,t) 

E i 

t n (x,e) 

t n ( X 'U x ) 

j x (x,t) 
A + 


Coulomb single scattering cross 
section 

Number transmission factor 

Cosine interval 

Current 

Initial or incident energy 
Function representing energy 
spectrum 

Function representing angular 

distribution 

Transmitted current 


INTRODUCTION 

When traversing even a thin layer of 
matter, electrons engage in numerous col- 
lisions which produce in most cases small 
energy losses and deflections. In addi- 
tion, they may undergo a relatively small 
number of catastrophic collisions which 
cause them to lose an appreciable fraction 
of their energy and to scatter through a 
large angle. The combined effect of all 
collisions is a complex transport process 
which requires an elaborate theory for de- 
scription. 

Monte Carlo methods of calculation 
have been applied to electron shielding 
calculations and in fact have been consid- 
ered up to now the most accurate avail- 
able, even though significant limitations 
are recognized. As with other Monte Carlo 
based calculations, the answers obtained 
can only be known with an accuracy which 
is governed by the statistical uncertainty 
inherent in the stochastic nature of the 
Monte Carlo method. Since individual 
electron slowing down case histories can 
be extremely complex, enormous amounts of 
computer time may be required to generate 
a statistically representative number of 
individual histories. This is the case 
even with the application of special tech- 
niques adapting Monte Carlo to the re- 
quirements of electron slowing down. 

Another problem with the Monte Carlo ap- 
proach is that it seems suitable for han- 
dling only steady state or time indepen- 
dent phenomena: at least the present au- 
thors are unaware of any charged particle 
Monte Carlo calculations structured to 
take into account time dependence# 

A computer-based method of treating 
the motion of charged or neutral particles 
which overcomes these difficulties has 
been introduced by Tavel and Zucker Cl, 2, 
3,4]. This technique, referred to as the 
"Phase Space Time Evolution" (PSTE) method, 
has been successful in the several neutron 


transport problems it has thus far been ap- 
plied to [4]. The present paper will deal 
with the first application of the PSTE 
method to a charged particle problem, 
namely, electron slowing down and shielding 
in a semi-infinite medium. 

The classical analytical approach to 
transport calculations usually requires the 
solving of an integro-dif f erential equation 
subject to boundary conditions established 
by the problem of interest. The equation 
itself is a mathematical representation of 
the space-time evolution of a particle dis- 
tribution which has been derived from the 
application of continuity principles in 
phase space. Instead of utilizing an in- 
tegro— different ial equation and solving it 
by computer techniques, the PSTE method 
makes direct use of the computer and em- 
ploys its bookkeeping capacity to literally 
keep track of the time development of the 
phase space distribution. 

The PSTE method was extended signifi- 
cantly beyond the approach used for neutron 
problems detailed in reference 4 in order 
to investigate the transport of electrons. 
This involved the addition of energy as a 
third dimension of the phase space and the 
use of multiple interaction theories for 
particle transport. 

Defining the Three-Dimensional Digitized 
Phase Space 

The phase space within which the 
flight of electrons can be traced which is 
shown in Figure 1 has Cartesian axes repre- 
senting energy, E, position, X and cosine 
of the polar angle with respect to the X 
axis, u . The values of E range from the 
maximum energy considered in the problem at 
the origin to the lowest energy of inter- 
est. The values of X can vary between any 
desired one dimensional spatial limit. The 
direction cosine, |j , ranges from +1 to -1. 

The particle density in the phase 
space is stored as the number of particles 
per unit length in the element of phase 
space between coordinates and X^ with 

energy E K and direction cosine u X j* 

SOLVING A PROBLEM 

An initial particle distribution rep- 
resents the state of the phase space 
at time zero. The time iterative scheme 
used in the calculation then traces the 
movement of each element of phase space 
for one time interval At. As illustrated 
in Figure 2 during this time interval the 
element of phase space is first relocated 
in the E-X plane. 
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Figure 1 An illuetretion of the three dimensional phaae epace, 


Mx 



In the calculations presented in this 
paper continuous slowing down theory as 
formulated by Rohrlich and Carlson [5], 
has been used to determine the energy loss 
of an element during At. Therefore, if 
the particles represented by an element of 
phase space are assumed to be traveling at 
an initial speed, v„, corresponding to 
their initial kinetic energy, e^, after a 
small time interval. At, their final ki- 
netic energy, Ej^', is approximately given 
by 

E k' = Ek " l dE k/ ds l v k At (1) 


where v 

v k = c y- - CmocV^+m^ 2 )] 2 ^ (2) 

and |dE k /ds| is the stopping power given 
in reference 5. By using the fact that 
the radiative contribution to slowing down 
is approximately given by [6] 

-dE k /dS (ZE k /800), ZE k «800 

where -dE k /dS is the energy loss per unit 
pathlength due to nonradiative collisions, 
the energy loss from radiative collisions 
has been included in eq. (1) by multiply- 
ing |dE k /ds| by the factor [l+ZE k /800]. 

It should be mentioned here that the PSTE 
method does not require the use of con- 
tinuous slowing down* Straggling (energy 
loss fluctuations) can also be handled 
within the framework of this method. It 
was intended, however, in the present ap- 
plication to first see how well the method 
performed without this complication, which 
has been left for future work. 

The technique employed to take into 
account the functional dependence on en- 
ergy of dE/dS and v in eq. 1 divides the 
time interval At into many smaller inter- 
vals during which it is assumed that dE/ 
dS and v remain constant. After each time 
increment the resulting final energy, E^', 
calculated is used to establish the values 
of dE/dS and v to be used for the next 
smaller time step. The total pathlength 
traveled during At and the final energy 
after At are obtained by* adding the con- 
tributions of all the smaller time steps. 

In addition to providing the total 
pathlength traveled and the final energy, 
the iterative procedure outlined for eq. 

(1) also generates afi average speed v k for 
the time interval At. Once this value is 
available the displacement of the phase 
space element can be found through use of 
the relation 

AS x = *k » i x j At (3) 


where Uxj is the initial direction cosine 
of the element and As x is its displacement 
after At. Once the X displacement has been 
determined the new position of the element 
X' is given by 

X* = X + AS X .. (4) 
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An important assumption in the determina- 
tion of the X displacement is that p x . does 
not change appreciably during At. This 
approximation is very good for small values 
of At. A constant Ux . during At is assumed 
mainly to avoid excessive computer time 
which contributes little to the accuracy of 
the electron PSTE calculation. 

After an element of phase space has 
been moved in the E-X plane, it must be 
apportioned to the established grid lines 
(digitized) . A moments weighting technique 
is used by the application presented in 
this paper to carry out the apportionment 
in energy. 

The procedure followed first estab- 
lishes the difference between the final 
energy and the end points of the energy 
grid interval into which the phase space 
element being operated upon falls. Thus, 
if some final energy of an element, E 7 , 
falls between E^ 7 -l and E^ 7 , these dif- 
ferences are E^'-l - E 7 and E 7 - E^ 7 , re- 
spectively. Once these values are found 
they are divided by AE , the energy inter- 
val between grid lines, to establish the 
energy apportioning fractions (moments). 


and 


v 


E k 


'-1 

-SIT 


e' 


f k ' -1 = 



(5) 

( 6 ) 


A similar procedure is involved in 
digitizing the final position coordinate. 
Suppose that the original element inter- 
sects the grid line ' with its end points 
falling on some and X 2 >X i / respec- 

tively. Then the distance X^i-X^ and 
X. , divided by the increment of distance 
between grid lines, AX, are the position 
apportioning fractions, 

F , = x i ’ ~ X 1 (7) 

AX 

and 

F ., +1 = X 2 ~ X i‘ . ( 8 ) 

AX 


The four fractions, and F^ ■ 

representing energy apportionment and F^« 
and + representing spatial distribu- 
tion, determine the digitization or appor- 
tionment in the E-X plane. Therefore, the 

number of particles per unit length con- 
tributed by the transported element having 
an initial particle density of p.-j^ j ^ 
(neglecting for the moment angular depen- 
dence) 


at 

E k'-1' 

between and X^ 1 

is 

p k'-l. j,i’ 

= Pk, j,i »k'-l F i ,]; 

at 


between X^ , and X^, + ^ 

is 


= p k, j,i tF k‘-l F i ’+1] ; 

at 

E k’ ' 

between X^._^ and X^, 

is 

p k\ j.'i* 

= Pk,j,i [ p k ,F i ,]f and 

at 

E k‘' 

between X^ » and X^, + ^ 

is 

p k\ 

= p k, j, i ^ F k' ,p i *+l^ * 


The only operation remaining to com- 
plete the movement of the original element 
for At is its distribution in u x - To per- 
form this operation the Goudsmit- 
Saunderson distribution [7,8] for multiple 
scattering was modified to provide the 
angular distribution of scattering refer- 
enced to the X direction independent of 
azimuthal angle, 

P(u x lui> = [jf V s ’> ds '] 

(9) 

where 

G^(s) = 2ttN^ ct(9 s .S) |l-P^(cos 9 s ) 

sin 0 d 0 ( 10 ) 

s s 

and u and 4 . are the cosines of the polar 
angles of scattering and incidence, re- 
spectively, referenced to the X direction. 
To utilize the modified angular distribu- 
tion in the PSTE calculation, the proba- 
bility density given in eq. (9) is eval- 
uated in the following way. First the in- 
terval - 1 % ^1 is divided into subinter- 
vals in sucJi a manner that a grid line 
lies in the middle of each interval. The 
probability, therefore, of obtaining some 
U x , in the interval is 9 iven b Y 

u a 

J = h ( Ua -u b ) 

Mb 
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Eq. (11) is applied to all subintervals. 

The probabilities generated in turn serve 
as weighting functions for the apportion- 
ment of the original transported element 
in u . 

Once every element in the phase space 
has been operated upon as outlined above, 
the stored results are tested against a set 
of terminating criteria for the problem 
being investigated. If these criteria are 
satisfied the calculation is completed. If 
not, another iteration is required. The 
final result of the PSTE calculation is a 
complete record of the time evolution of a 
particle density distribution given as a 
function of energy, direction and position. 


Comparison of PSTE Calculations with Elec- 
tron Monte Carlo and Experimental Results 


To assess the validity of the results 
provided by the PSTE method as applied to 
electron transport, a comparison was made 
with Monte Carlo calculations and experi- 
mental results for a shielding problem 
which considered a 1 MeV pulse of electrons 
normally incident on a semi-infinite slab 
of aluminum. Comparisons were made on the 
basis of the different forms of published 
Monte Carlo and experimental results for 
this problem. 

The first form of result compared is 
the number transmission factor. This 
factor, T n (X), which provides the fraction 
of incident electrons transmitted past X, 
is defined by 


1 t 


* AEi U 


k HI 

J o J o J o 


J(X,E,u x ,t')dt , d Ux dE 


( 12 ) 


where the current J(X,E, Ux ,t) is given by 


J(X,E /k i x ,t) = p(X # E,n x ,t)v(E)u x , (13) 


E^ is the incident energy of the electrons 
and t is the time required for the entire 
pulse to essentially pass through the me- 
dium. Comparison of the PSTE generated 
number transmission curve with several 
Monte Carlo and experimental results is 
given in Figure 3. Penetration (X) in this 


case is represented as the fraction of the 
mean range traversed. The Monte Carlo re- 
sults include calculations by Berger [9,10, 
11] which use continuous slowing down in 
one case and straggling (energy loss fluc- 
tuation) in the other. Also shown are 
Monte Carlo calculations made by Perkins 
[12] which include the effects of strag- 
gling. The experimental points depicted 
include the results of experiments con- 
ducted by Rester [13] for 1 MeV electrons. 
Since the number transmission curve plotted 
as a function of the mean range traversed 
is approximately independent of initial en- 
ergy in the neighborhood of 1 MeV, the av- 
erage results of Agu et al. [14] for ex- 
periments conducted at. energies below 1 
MeV have also been used for comparison pur- 
poses. The agreement exhibited in Figure 3 
between the PSTE and the Monte Carlo and 
experimental results is good. The PSTE 
curve, as should be expected, falls below 
the Monte Carlo and experimental results 
which considered straggling effects since, 
as mentioned before, this PSTE calculation 
utilized the continuous slowing down ap- 
proximation for energy degradation. Al- 
though the experimental points of Rester 
lie somewhat above the PSTE curve, it 
should be noted that the PSTE values fall 
within the range of experimental uncer- 
tainty associated with these results. 



Figure 3 Comparison of a PSTE number transmission curve with Monte Carlo and experi- 
mental results for a 1 MeV pulse of electrons normally incident on a semi- 
infinite slad of aluminum. 
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The energy spectrum of transmitted 
electrons produced by the PSTE calculation 
for the problem under consideration is a 
critical indication of agreement with ex- 
isting calculations and experimental re- 
sults, The energy spectrum T n (X # E) is re- 
lated to the PSTE generated current by 

1 t 

v x - E ’ --dzff 

x J o J o 

J(X,E, Ux ,t')dfd Ux • (14) 

The shapes of PSTE energy spectra for 
three different penetrations in terms of 
fractions of the mean range are compared to 
Monte Carlo and experimental results docu- 
mented by Rester [13] in Figures 4-6. 

The Monte Carlo calculations were performed 
by Berger. Assuming that the total number 
of particles present at each penetration 
investigated is approximately the same for 
the PSTE, Monte Carlo and experimental ex- 
amples, the histograms in Figures 4-6 
have been drawn on a scale relative to the 
maximum of each distribution. In this way 
it is possible to separate an examination 
of the energy spectrum from other consid- 
erations. In other words, it is only in- 
tended here to investigate the relative 



Figure 5 


Comparison of the PSTE transmitted energy apectrum with Monte Carlo and ex- 
perimental results at approximately -4 mean range for a 1 MeV pulse of elec- 
trons normally incident on a semi-infinite elab of aluminum. 



Figure 4 Comparison of ths PSTE trsnsmittsd snsrgy spectrum with Monte Carlo and ex- 
perimental results at approxlmstsly .2 maan rang# for a 1 MeV pulee of elec- 
trons normally incidsnt on a ssmi-inf inite slab of aluminum. 



Figurs 6 comparison Of the PSTE transmitted energy spectrum with Monts Carlo and ex- 
perimental results at approximately .6 mean range for a 1 MeV pulae of elec- 
trons normally Incident on a semi-infinite slab of aluminum. 
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shapes of the different energy spectra and 
not the total number of particles at a 
specific • penetration (this has been done in 
the consideration of the number transmis- 
sion factor). The good general agreement 
between Monte Carlo and experimental re- 
sults and the PSTE calculation is well il- 
lustrated in Figures 4-6. For each pene- 
tration the PSTE result agrees very well 
with the Monte Carlo histogram. The agree- 
ment is not as good when the PSTE results 
are compared to the experimental values. 

The largest discrepancy appears at .4 of 
the mean range where the disparity between 
the Monte Carlo and experimental results 
is the greatest. 



ANGLE, DEGREES 


Figure 7 Comparison of the PSTE transmitted angular distribution with Monte Carlo and 
experimental result, at approximately .2 mean range for a 1 MeV pule® of 
electron, normally incident on a aemi-inf inite .lab of aluminum. 



ANGLE, DEGREES 


Figure 8 Comparison of the PSTE transmitted angular distribution with Monte Carlo and 
experimental results at approximately .4 mean range for a 1 MeV pulae of 
electrons normally incident on a semi-inf inite slab of aluminum. 



ANGLE, DEGREES 

Figure 9 Comparison of the PSTE transmitted angular distribution with Monte Carlo and 
experimental results at approximately .6 mean range for a 1 MeV pulae of 
electrons normally incident on a semi-infinite slab of aluminum. 


The angular distribution of transmit- 
ted electrons produced by the PSTE calcula- 
tion has also been compared to Monte Carlo 
and experimental results. The angular dis- 
tribution T n (X, |a x ) is defined by 



0 0 

J(X,E #Ux ,t')dt'dE (15) 

where 0 £ The shapes of the angular 

distribution curves for the same three 
penetrations used in the energy spectrum 
investigation are compared in Figures 7 - 
9. Again the curves are drawn on a scale 
relative to the maximum of each distribu- 
tion for the same reasons outlined in the 
discussion of the energy spectrum. The 
results used for comparison are Monte Carlo 
calculations made by Berger and experiments 
conducted by Rester, the data for both 
being taken from publications authored by 
Rester [13,15]. In general there is ac- 
ceptable agreement between the PSTE re- 
sults and the Monte Carlo and experimental 
points. For all three penetrations the 
PSTE curve agrees very well with the ex- 
perimental values. The largest discrep- 
ancy between the PSTE and Monte Carlo re- 
sults occurs at the smallest and largest 
penetrations considered. Interestingly 
the greatest difference between the Monte 
Carlo and experimental values also occurs 
at these depths. 
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Figure 10 Transmitted current at 0.22 mean range for a 1 MeV pula* of electrons nor- 
mally incident on a semi-infinite slab of aluminum resulting in an incident 
current of about 3.2 X 10 12 electrons per second. 



Figure 11 Transmitted current at 0.4 mean range for a 1 MeV pulae of electrons nor- 
mally incident on a aemi-inf inite slab of aluminum resulting in an incidsnt 
current of about 3.2 x 10 12 electrons per second. 



Figure 12 


Transmitted current at 0.62 mean range for a 1 MeV pulse of electrons nor- 
mally incident on a semi-infinite slab of aluminum resulting in an incident 
current of about 3.2 x 10 12 electrons per second. 


Time Dependent Results 

One of the more significant contribu- 
tions of the PSTE method is the time de- 
pendent results it provides. In Figures 
10 12 plots of the time dependent trans 1 

mitted current J v (X,t) given by 
x + 

E 

V x ' c) 

0 

J(X,E,M x ,t)d Ux dE (16) 

for the problem in question are presented. 
The penetrations considered are the same 
as those examined for the energy spec- 
trum and angular distribution. At the 
small penetration* (.22 mean range) 
the transmitted current reaches a max- 
mum very quickly and then decreases 
rapidly with time. The intermediate pene- 
tration (.4 mean range) results exhibit 
much of the same behavior but the current 
drops off at a slower rate as time in- 
creases. At the largest penetration (.62 
mean range) there is a more gradual build- 
up to the maximum transmitted current 
achieved and a less rapid decrease with 
time than at the other penetrations. This 
is probably because at this depth the elec- 
trons are clif fusing through medium and are 
characterized by relatively broad energy 
spectra and angular distributions. 


/ 
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Advantages of the PSTE Method 


Another way of looking at the time 
dependent results provided by the PSTE 
method is illustrated in Figure 13. Here 
the current as a function of position with 
time as a parameter has been traced by the 
computer. The time interval between plot- 
ted iterations is approximately 3.2 x 
1CT 13 seconds. Essentially this figure 
represents a picture of the time evolution 
of the transmitted current in the aluminum 
slab. At first the current builds up rap- 
idly and declines rapidly as the electrons 
enter the medium. There is also a signif- 
icant reduction in the maximum transmitted 
current from one iteration to the other. 

At larger penetrations or at a later time, 
however, the curves representing the dis- 
tribution of the transmitted current as a 
function of penetration become broader. 

This is understandable due to the fact that 
at later points in time the electrons are 
diffusing through the medium and have broad 
energy spectra and angular distributions. 
This parallels the reasoning used to de- 
scribe the shape of the transmitted cur- 
rent curve as a function of time at the 
largest penetration examined (.62 mean 
range. Figure 12) . 


One of the most attractive features 
of the PSTE method is the comparative J.y 
small amount of computer time it requires. 
For example, the PSTE calculations dis- 
cussed in this paper required approximately 
ten minutes of computer time on the CDC 

♦ 

6600 computer at the Brookhaven National 
Laboratory. This compares to a Monte 
Carlo calculation for the same problem 
which requires on the order of an hour or 
more of computer time. 

Another advantage of the PSTE method 
is that it provides deterministic results. 
This can be contrasted to the results of a 
Monte Carlo calculation which have a sta- 
tiscal uncertainty associated with them 
dependent on the number of histories 
sampled. 

Finally, the PSTE method generates a 
complete record of the time evolution of a 
particle density distribution. As stated 
previously, the authors are unaware of any 
successful Monte Carlo attempts at provid- 
ing time dependent results for charged 
particle problems. It is the authors' 
opinion that the PSTE method is ideally 
suited for this application. 



Figura 13 Ti». .volution of transmitted curr.nt for a 1 HeV pulse of alactrona nor- 
mally incident on a semi-inf inita slab of ^ltninta* raaulting In an incid.nt 
currant of about 3.2 x 10* 2 alactrona par aacond. 
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